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The bond order wave (BOW) phase of the extended Hubbard model (EHM) in one dimension (ID) 
is characterized at intermediate correlation U — At hy exact treatment of A^-site systems. Linear 
coupling to lattice (Peierls) phonons and molecular (Holstein) vibrations are treated in the adiabatic 
approximation. The molar magnetic susceptibility Xm{T) is obtained directly up to = 10. The 
goal is to find the consequences of a doubly degenerate ground state (gs) and finite magnetic gap 
Em in a regular array. Degenerate gs with broken inversion symmetry are constructed for finite N 
for a range of V near the charge density wave (CDW) boundary a,t V ~ 2.18f where Em ~ 0.5t 
is large. The electronic amplitude BiV) of the BOW in the regular array is shown to mimic a 
tight-binding band with small effective dimerization 5e//. Electronic spin and charge solitons are 
elementary excitations of the BOW phase and also resemble topological solitons with small (5e//. 
Strong infrared intensity of coupled molecular vibrations in dimerized ID systems is shown to extend 
to the regular BOW phase, while its temperature dependence is related to spin solitons. The Peierls 
instability to dimerization has novel aspects for degenerate gs and substantial Em that suppresses 
thermal excitations. Finite Em implies exponentially small Xm{T) at low temperature followed by an 
almost linear increase with T. The EHM with U = At is representative of intermediate correlations 
in quasi-lD systems such as conjugated polymers or organic ion-radical and charge-transfer salts. 
The vibronic and thermal properties of correlated models with BOW phases are needed to identify 
possible physical realizations. 
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I. INTRODUCTION 



Nakamura identified the bond order wave (BOW) 
phase [H of the one-dimensional (ID) half-filled ex- 
tended Hubbard model (EHM, Eq. [J). The key 
features are broken inversion symmetry, doubly degen- 
erate ground state (gs) and finite magnetic gap Em in 
a regular (equally spaced) array. Competition among 
electron delocalization t, on-site repulsion U > and 
nearest-neighbor repulsion > stabilizes the BOW 
phase over a narrow range whose boundaries motivated 
subsequent studies 0-11]. The BOW phase has V « U/2 
and substantial t. The quantum transition to the charge 
density wave (CDW) phase at large V is first order 
i, i] for f7 > f/* w 7t, continuous for U < U* . The 
BOW/CDW boundary is at Vc{U) for U < U* , while 
the boundary to the spin-fluid phase with = is 
at Vs{U) < Vc{U). Other half-filled ID Hubbard mod- 
els with spin- independent interactions have narrow BOW 
phases between the spin- fluid and CDW phases p^ . 

Theoretical studies of the EHM have focused on the 
quantum phase diagram of an extended array without 
reference to possible physical realizations. Even for mod- 
els, however, ID instabilities and finite temperature must 
be addressed. The principal goal of this paper is to char- 
acterize the BOW phase of the EHM at U = At, a typ- 
ical choice for intermediate correlation, with special at- 
tention to the consequences of gs degeneracy and finite 



Em- Coupling to Peierls or Holstein phonons requires gs 
derivatives in the adiabatic (Born-Oppenheimcr) approx- 
imation, and finite temperature properties are needed 
to assess physical realizations. It is clearly desirable to 
understand BOW phases prior to specific applications. 
Broadly similar properties are expected for other U or 
other ID models with spin-independent interactions. 

Quasi- ID organic molecular crystals and conjugated 
polymers have strong electron-phonon (e-ph) coupling 
and clectron-molecular-vibration (e-mv) coupling in ad- 
dition to intermediate correlation. The Su-Schrieffer- 
Heeger (SSH) model of polyacetylene has linear e-ph cou- 
plingand topological solitons as elementary excitations 
[111 . KL^ . The Peierls instability is driven by e-ph cou- 
pling and has spectacular e-mv consequences in infrared 
spectra when inversion symmetry is broken. Conjugated 
polymers and organic ion-radical salts have been a play- 
ground for ID Hubbard models p^ - [T7| with e-ph and 
e-mv coupling, variable electron or hole filling, degener- 
ate or nondegenerate gs, and either segregated or mixed 
stacks. The bandwidth is At ~ 10 eV in polymers and 
4t « 1 eV in 7r-stacks, with comparable intermediate t/U . 
The Peierls instability of Hubbard-type models is a rich 
separate topic 18] . Wc arc not aware of work on either e- 
ph or e-mv coupling in the BOW phase. As shown below, 
finite Em and degenerate gs lead to electronic solitons in 
a rigid regular array that nevertheless resemble SSH soli- 
tons. With suitable modification, extensive SSH analysis 
p^ . [l^ can be applied to models with a BOW phase. 
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We recently proposed that a BOW phase is reahzed 
in Rb-TCNQ(II), the second polymorph of a tetra- 
cyanoquinodimethane salt [l^, [l^l- The evidence is a 
100 K crystal structure (Pi) with regular TCNQ" stacks 
at inversion centers, negligible spin susceptibility below 
140 K that indicates a large Em, and infrared spectra 
that demonstrates broken electronic inversion symmetry. 
Broken Ci symmetry and finite Em in a regular ID 
array are precisely the signatures of a BOW phase [l|. 
Large E^ indicates proximity to the CDW boundary, 
and we will so choose V in the EHM a.t U = At. We have 
Rb-TCNQ(II) and alkali-TCNQs in mind, but do not 
model them explicitly beyond invoking solitons for the 
temperature dependence of infrared spectra. In addition 
to values of microscopic parameters, interactions be- 
tween chains must be addressed in actual models along 
with the Coulomb interactions and transfer integrals for 
different stacking motifs. We consider electronic proper- 
tics of the EHM at U = At with linear coupling to Peierls 
and Holstein phonons. Vibrational degrees of freedom 
are introduced as needed in the adiabatic approximation. 
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The EHM describes electronic degrees of freedom in a 
regular ID array 

N 

Hei = ^ -t{al,^ap+i^„ + h.c) 

p—l,a 
N 

+ ^{Unpinp-l)/2 + Vnpnp+i) (1) 
p=i 

where i = 1 is the unit of energy, h.c. is the hermitian 
conjugate, a| g.(apo-) creates (annihilates) an electron 
with spin a at site p, and Up is the number operator. 
H conserves total spin. The gs is a singlet {S = 0) for 
U, V > and even N. The half-filled band with N 
electrons and N sites has electron-hole (e-h) symmetry 
J = ±1 and inversion symmetry Ci at sites that we 
label as a = ±1. The correlated many-electron basis 
increases as « 4^ at large N. Valence bond (VB) 
methods [2l|, at present yield exact results up to 
= 17 for low-energy states and the full spectrum up 
to N = 10. Nakamura identified [l| the BOW phase 
using field theory, symmetry arguments and numerical 
results up to N = 12. 

To introduce the EHM phase diagram at K, we 
define three threshold excitations from the singlet gs: 
the magnetic gap Em to the lowest triplet, the gap Ej 
to the lowest singlet with opposite J, and the gap E^ 
to the lowest singlet with opposite Q. Increasing V at 
constant U, t drives the system from a spin-fluid phase 
to a CDW phase. Fig. [1] shows the evolution of Em, 
Ej and E^ with V for U = At, N = 12. The BOW 
phase spans Vg < V < Vc, where Vs{N) and Vc{N) are 
defined by the excited-state crossovers Ea{N) = Em{N) 
and EaiN) = Ej{N), respectively. Table |T] lists the 



FIG. 1: Excitation thresholds and crossovers of the 12-site 
extended Hubbard model, Eq. [T] with periodic boundary con- 
ditions. Em, Ea and Ej are energy gaps to the lowest triplet 
and the lowest singlets with opposite d and e-h symmetry, 
respectively. Vs is defined by Ea = Em, Vi by Ecr ~ and Vc 
by E, = Ej. 



TABLE I: Excitation thresholds and crossovers of the ex- 
tended Hubbard model, Eq. [T] with A'' sites, U — A, and 
tiN = ±1 for iV = An, An + 2. 
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remarkably weak size dependence of Vs{N) and Vc{N) 
for N = An with periodic boundary conditions (PBC, 
tiN = 1) and N — An + 2 with antiperiodic boundary 
conditions {tiN = ^1)- Similarly, the Vg boundary 
of a frustrated spin chain has been found |.23j using 
Em{N) = E^{N) up to iV = 24. Multiple methods yield 
the boundary Vc{N) of the CDW phase 0, The 
points V = Vi {N) in Table [J or Fig. [T] correspond to 
Ecr — 0. They mark a gs degeneracy that is central to 
our discussion, where broken symmetry gs are readily 
constructed. 

The paper is organized as follows. Broken inversion 
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symmetry and elementary excitations are treated exactly 
in Section II for finite N. Electronic solitons, both spin 
and charge, are found in regular chains with open bound- 
ary conditions (OBC, tiN — 0) and compared to SSH 
solitons. Section III deals with linear coupling to molecu- 
lar (Holstein) and lattice (Peierls) vibrations. The Berry 
phase formulation of polarization is applied to the in- 
frared activity of molecular vibrations when Ci symme- 
try is broken. The Peierls instability of the BOW phase 
is contrasted to the SSH model within the limitations of 
an adiabatic approximation. The magnetic gap and 
spin susceptibility xm of the BOW phase are obtained 
in Section IV for large Em close to the CDW instabil- 
ity. Large Em reduces the thermal population of excited 
states and opens a new regime in which spin solitons 
govern xm ■ The discussion in Section V summarizes the 
consequences of broken symmetry and finite Em such as 
the temperature dependence of the infrared intensity or 
of Xm- We briefly mention extensions to BOW phases of 
related models. 



II. BROKEN SYMMETRY AND ELEMENTARY 
EXCITATIONS 

We discuss the EHM, Eq. [T]with U = 4, t = 1, using 
exact results for finite N with PBC(tiAr = 1) for N = An 
and tiTv = —1 for N — in + 2. The gs kinetic energy is 
given by the bond orders p„ of successive sites 



that is not available for Monte Carlo, which yields 
the gs energy, or density matrix renormalization group 
(DMRG) calculations. DMRG with OBC breaks d 
symmetry for even N and returns a nondegenerate gs. 
We are not aware of DMRG with PBC that conserves Ci 
and gives the energy and gs in both the a — ±1 sectors. 
It is easier to find an energy crossover at Vi{N) where 
Ea- — for finite N than to demonstrate a degenerate gs 
over the interval Vg < V < Vc in the extended system. 

On the other hand, we need the gs degeneracy at U = 
At beyond the single point Vi{N). To do so, we add the 
following perturbation to Eq. [1] 

H{J2)^J2j2Sp.Sp+2. (5) 

p 

H{J2) acts on second neighbors with one electron each. 
The upper panel of Fig. [5]shows how Fi( J2) scans E„ 
across the BOW phase at C/ = At, lowering Vi for antifer- 
romagnetic J2 > and raising it for J2 < 0. The interval 
J2 = ±0.15 in Fig. [5] is sufficient to enforce Ea{J2) — 
between V/t = 2.0 and 2.15 for the EHM with = 12 
and U/t — A. Strictly degenerate gs in the cr = ±1 sec- 
tors will be essential for electron vibrational coupling in 
Section HI. 



2pn = {ipo\^{al,^an+i,a + h.c)\ipQ) 



(2) 



Broken Ci symmetry in the BOW phase leads to p2n 7^ 
P2n-i , while broken e-h symmetry in the CDW phase 
leads to different electron count n2p n2p-i in the even 
and odd sublattice. At constant U and t, the order pa- 
rameter B{V) of the BOW phase is 



BiV)^\p2n{V)-p2n-liV)\. 



(3) 



B{V) is large between sites 2n, 2n — 1 in one broken- 
symmetry gs and between 2n, 2n + 1 in the other. 

Since both Em and E^ vanish rigorously in the spin- 
fluid phase with V < Vs, finite gaps in Fig. [1] are due 
to finite N. The dashed hne Em{V) — E„{y) is an ap- 
proximation for opening the magnetic gap. Similarly, the 
dashed line Ej{y) — E^^y) approximates the closing of 
the e-h gap at Vc- Finite N limits Ef, = to points 
V\{N) in Fig. [T]and Table ID At gs crossovers, we con- 
struct broken-symmetry states 



l^±(Vi)) = (|^.=i)±|^.=-i))/\/2 



(4) 



and compute their bond orders in Eq. [21 Degenerate 
gs at finite N provide direct access to BOW systems 
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FIG. 2: (top panel) Energy of the lowest singlet with 
opposite Ci symmetry in a 12-site EHM with U = At and PBC 
as a function of V in Eq . [l]and the indicated J2 values in Eq. 
[5l (bottom panel) Order parameter -B(Vi, J2) = |p+ — P-| in 
Eq. [3] at S,^ = for iV = 12 (closed symbols) and B{V) for 
iV = 8, 12, 16 for J2 = (open symbols). 
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The gs in Eq. |4] and the order parameter B in Eq. 
[3] hold for Vi{J2) where E„ — 0. The lower panel of 
Fig. [2] compares B{Vi{J2)) at the indicated J2 values 
for TV = 12 with B{V) based on J2 = 0, when \tp±) 
in Eq. |4] are linear combinations of functions that 
are not quite degenerate. B{V,N) depends weakly 
on and N in this interval. Of course, H{J2) 
also perturbs other properties, at least slightly, and 
becomes a strong perturbation for V < 2t. We are inter- 
ested in V close to Vc, where is small and E^ is large. 

Bond orders p±{V,N) and B{V,N) are well approxi- 
mated near V ~ Vi{N) without invoking J2, and results 
in this Section are based on J2 =0. For example, U — At, 
V = 2mt, = 16 returns p+ = 0.752, p_ = 0.398, 
B — 0.354. As expected [l^] near the metallic point Vc, 
the average p = 0.574 is within 10% of 2/tt, the value 
for a tight binding (Hiickel) band of free electrons with 
U — V = in Eq. [TJ The t term dominates near Vc 
where U and V almost cancel. 

The SSH model [III El is a tight-binding band with 
linear e-ph coupling a = {dt/du)Q to the Peierls phonon, 
the optical mode of the ID chain, and an adiabatic ap- 
proximation for a harmonic lattice. Its gs is dimerized, 
with tn = — (1 — (5(— 1)") along the stack. Bond orders 
of the infinite chain are readily found analytically, with 
p+{S) = 0.806, p-{6) = 0.452 and B{6) = p+ ~ p_ = 
0.354 at (5 = 0.10 that matches B{V) of the finite EHM 
above. The origin of broken Ci symmetry is quite dif- 
ferent: e-ph coupling in SSH, correlation U,V in EHM. 
Moreover, B{V) opens at Vg « 1.86i and is almost con- 
stant for V > 2.0t before vanishing abruptly at Vc ~ 2.18t 
while B{6) is monotonic in 6. The electronic gs are nev- 
ertheless similar and previous discussions of topological 
solitons or domain walls between regions with opposite 
B can be applied to BOW phases [ll [ll [H . 
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FIG. 3: Bond orders pr in Eq. [2] of a spin soliton (close 
symbols) or charge soliton (open symbols) in an A'^-site EHM 
with tiN = 0, U = 4t andV ^ 2.05f in Eq[Tl 



U = V — and odd N can be read as a Hiickel model 
of an alternant hydrocarbon with a nonbonding orbital 
at e = that is empty in the cation, singly occupied in 
the radical and doubly occupied in the anion. Since the 
e = orbital has nodes at every other site for any N, it 
does not contribute to bond orders that are consequently 
equal for neutral and charge solitons in the SSH model. 
The correlated model has equal Pr for charge solitons 
that differ slightly from the spin soliton. In every case, 
bond orders are p+ at the ends and reverse smoothly in 
between. 



We start with electronic domain walls and return later 
to the Peierls instability. Domain walls are modeled 
as in SSH with odd N in Eq. [T] and tiN = 0. We 
retain PBC for V to minimize end effects. Ne = N 
electrons correspond to a neutral soliton with S — 1/2, 
while iV ± 1 electrons are charge solitons with 5* = 0. 
Spin-charge relations are reversed just as in SSH. We 
label sites from r = at the center to r = zL{N — l)/2 
at the ends. Positive and negative solitons have equal 
bond orders by e-h symmetry. We could solve = 17 
exactly using a symmetry- adapted valence bond basis 
to compute p„ in Eq. [2l As shown in Fig. [3] for 
= 15 and 17, the bond orders of spin and charge 
solitons are symmetric about the center. They are 
slightly larger at the chain ends than p_ ~ 0.398, 
p+ = 0.752, the values of the BOW phase at V ^ 2.05*. 
This end effect for tiN = also appears for even A^ 
and alternating bond orders with largest p+ at either end. 

As context for these results, we note that Eq. [T]with 



Next we compute spin densities pr — 2{S^.) for the 
radical and charge densities qr = 1 — 11,. for ions. 
The upper panel of Fig. 2] has pr for A^ = 15 and 
17 while the lower panel has qr for the cation; the 
anion charges are —qr by e-h symmetry. The spin 
or charge density vanishes in alternant hydrocarbons 
where the nonbonding orbital has nodes, at odd r for 
A^ — 4?! -I- 1 and even r for A^ = 4n — 1. As seen in 
Fig. m correlation 26] generates small negative pr at 
these sites and also small qr of opposite sign. The spin 
or charge density is large in the middle and decreases at 
the ends. Electronic solitons in the BOW phase differ 
in this respect from a regular tight-binding band, which 
has equal pr or qr at every other site for any odd A^. 
To mimic the BOW results in Fig. |4]with SSH solitons, 
finite 5 « 0.05 is needed for pr and 5 « 0.08 for qr 
since large 5 gives faster decrease. Electronic solitons in 
Figs. [3] and 13] connect regions with ±B in a regular chain. 

Domain walls or topological solitons are the elemen- 
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FIG. 4: (a) Soliton spin density 2(5^) and (b) soliton charge 
density q,- — 1 — rir of an A''-site EHM with tiN = 0, C/ = 4t 
and V = 2.05t in Eq. [T] 



tary excitations of the BOW phase. They resemble SSH 
sohtons in an appropriately dimerized lattice. The BOW 
phase has different energy for creating a pair of spin or 
charge solitons. We take Eo{N,Ne) as the gs energy of 
Eq. [l]for even N, Ne electrons and tiN = 0. The forma- 
tion energy of a pair of spin solitons (S — 1/2) is 



2WsiN) = EQ{N+l,N + l) + EQiN-l,N-l) 

- 2EoiN,N). (6) 

In the limit of large N, parallel spins are the lowest triplet 
with 2Ws = Ejn, the magnetic gap for PBC. The results 
for = 8, 12 and 16 in Table [n] are for V = Vi{N) in 
Table HI where Ea- = 0, but the V dependence is weak. 
The formation energy for a pair of charge solitons is 

2Wc{N) ^Eo{N + l,N) + Eo{N-l,N)- 2Eo{N, N).{7) 

The cation, anion and neutral system are singlets with 
different charge. The natural comparison is to the charge 
gap, 



2Ec{N) ::^Ea{N,N + l) + Eo{N,N^l)- 2Ea{N, N).{8) 

For large N, Ec{N) is the energy of separated ion 
radicals. As seen in Table [III 2Wc approaches Ec from 
below, as expected since charge solitons are singlets 
while the ions in Eq are doublets. On the other hand, 
2Wc is larger than Ej because > favors adjacent 
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sites with n = (hole) and n = 2 (electron). 

Valence bond methods yield low-energy excitations in 
every symmetry sector with fixed S, J and a. Except 
for total wave vector k = or n, degeneracy in ±k is 
expected and found. Finite-size effects increase with en- 
ergy and it becomes progressively more difficult to ex- 
tract more than 3-4 states for large N. The entries in 
Table [ill are from a much larger set. An "effective" 
Seff ~ 0.05 — 0.10 is inferred from the order parameter 
B{V) or from spin or charge solitons in the EHM with 
[/ = 4i and 2.0 < V/t < Vc- While BOW-phase results 
are for a regular stack, they are naturally related to SSH 
results with 5 < 0.10 that is considerably smaller than 
(5 = 0.18 based on the optical gap of polvacetvlene [iTi . 



III. COUPLING TO HOLSTEIN AND PEIERLS 
PHONONS 

The SSH model invokes linear e-ph coupling, 
a = (di/c?u)o, to characterize the Peierls instability 
and elementary excitations of a half-filled tight-binding 
band. Linear coupling to molecular vibrations is the 
basis for interpreting polarized infrared spectra of 
TT-radical stacks. The operators aj,^., Up^a create, anni- 
hilate electrons with spin a in the lowest unoccupied 
molecular orbital of TCNQ, with equal energy A = 
and Up = 1 for a TCNQ~ stack. Charge fluctuations 
[131 modulate A and illustrate linear Holstein coupling 
gn to the n*^ totally symmetric (ts) molecular vibration. 
When Ci symmetry is broken, ts modes become strongly 
IR allowed by borrowing intensity from the optical 
charge-transfer excitation and they are polarized along 
the chain. Accordingly, the appearance of ts modes in 
polarized IR yields detailed information [28l - l30j about 
e-mv coupling constants gn and has been widely used 
to infer dimerization. Charge fluctuations break e-h 
symmetry, and strict degeneracy at Vi{N) is critical 
because A modulation is a small energy. 

The Berry-phase formulation [3l|, [s^l of polarization 
makes possible improved vibronic analysis of extended 
systems. It is directly applicable to quantum cell models 
[33| such as the EHM. The polarization P per unit charge 



6 



and unit length is a phase [32, 12: 
27r 

Pn = —Im{lnZN) 

Zn = {^o\exp{2TnM/N)\^jjo} 



N 



M 



(9) 
(10) 

(11) 



where |V'o) is the exact gs of an A^-site supercell and 
M is the conventional dipole operator for a regular 
array with unit spacing. |2^jy| = is a metallic point 
at which P is not defined jsij; it corresponds to Vc for 
U < U* and a continuous CDW transition [l^l- The 
EHM has real Z and P = by either C,; or e-h symmetry. 

The required generalization of Eq. [T] for Holstein cou- 
pling in the adiabatic approximation for molecular sites 
is [SJ 



H{A) = H,, 



N 

A^(-l)^ 
p=i 



(12) 



Finite A breaks e-h symmetry, but not Ci symmetry. 
The gs in the BOW phase is now |-0±(V,A)) with J2 
in Eq. [5] chosen to have = for A = 0. Strictly 
degenerate cr = ±1 ensures first-order correction in A. 
Finite A mixes in excited states such as Ej in Fig. [T] 
whose energy in the BOW phase decreases rapidly with 
increasing V. We again vary V aX U = A, t = 1. Aside 
from a multiplicative constant, the IR intensity goes as 

M 



I 



IRiV) - ( 



(13) 



Charge fluctuations give a finite derivative at A = 0. 
The IR intensity in Eq. [T3] is purely electronic. It can 
be partitioned among ts modes as discussed [28l - [30j in 
systems where dimerization breaks C'i symmetry. 

We compute the derivative in Eq. [T3] in the BOW 
phase at V = Vi(J2). The imaginary part of Z^ is 
initially proportional to A, as shown in the inset of Fig. 
El Iir{V) increases rapidly with V up to Vc, as expected 
for decreasing Ej, and vanishes abruptly in the CDW 
phase where Ci symmetry is restored and e-h symmetry 
is broken. The corresponding band result [Sj, [s^] for 
(9P((5, A)/9A)q increases as 1/S^ and diverges at the 
metallic point A = (5 = where, however, P is not 
defined. Since the real part of Zn{V) =OatV = VciN), 
the finite system also has divergent {dP{V, A)/dA)o 
at Vc{N) and the BOW phase again resembles a band 
with small 5, with two major differences. First, Iir(V) 




FIG. 5: Infrared intensity Iir{V), Eq. 1131 of molecular vi- 
brations due to A in Eq. [12] and broken electronic symmetry 
in the BOW phase of the EHM. The inset is the polarization 
P(A) in Eq. [5] Jjjj vanishes by symmetry for V > Vc or 

< Vs. 



increases with V up to Vc while B{V) is almost constant. 
Second, Seff for IR intensity decreases with increasing 
V up to Vc while deff for B{V) increases with V from 
Vs . The SSH model has a single band gap of ASt instead 
of the EHM's multiple threshold excitations in Table |T| 
and Fig. [Tj For example, the metallic point Vc has large 
E^ « OM. 

Domain walls introduce inversion centers that reduce 
IiR between regions with opposite B. More quantita- 
tively, we again consider systems with odd N, V = 2.05 
and ti]\f = 0. We break Ci symmetry with ±A at sites 
±r from the center and evaluate {dP/dAr)o. Since A in 
Eq. [TH applies all sites in Fig. [SJ we scale P{V,Ar) by 
N/2 for the soliton. The inset to Fig. |6| shows P{V, A^) 
at y = 2.05* for A^ = 0.05, with P = at r = and 
scaled P « 0.06 at the ends of a spin soliton. The same 
pattern is found for charge solitons (data not shown) 
with about 50% higher intensity. In either case, Iir{V) 
at the ends is an order of magnitude less than the Vi(J2) 
values for degenerate gs in Fig. [5] This finite-size effect 
can be traced to higher Ej in the radical. We did not 
solve TV = 17 since finite A doubles the dimensions of 
the many-electron basis. IR intensities are consistent 
with a soliton and an inversion center between regions 
with ±B{V). 

We consider next the Peierls instability. The SSH 
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0.08 




FIG. 6: Infrared intensity Iir{V), Eq. 1131 of molecular vi- 
brations due to A in Eq. [12] in 13 and 15-site spin solitons 
with V = 2.05t and tijv = in Eq. [1] The inset shows 
P{V, Ar) for Ar = 0.30 in Eq. |l]at site r, with r = at the 
center and r — ±{N — l)/2 at the ends. 



(5 for T > 0. 

Our discussion is limited to eT{S) in Eq. 1151 The gs 
of correlated models is unconditionally dimerized when 
the regular array has = and Xd{S) = —{d'^eo/dS'^) 
diverges at 6 = 0. Examples [13] include the ionic phase 
of organic charge transfer salts and the spin fluid phase of 
Hubbard models or Heisenberg spin chains. Dimerization 
is conditional at edXd(O) = 1 when Xrf(O) is finite in the 
neutral phase of CT salts or the CDW phase of the EHM. 

Degenerate gs is different, as seen for NeriS) in Fig. [7] 
for V = Vi{N) in Table U and inverse stiffness td = 0.12. 
Now eo(^) goes as -B\6\ - Xd(0)5V2. The cusp B[V) is 
due to degeneracy at (5 = while the energy gap i?3 in 
the singlet sector, listed in Table HIl ensures finite X£i(0). 
Minimization of et yields 



8,^ - ±e<iB(F)/(l - edXd(O)) (16) 

with beq > for (5 > and (5eq < for (5 < 0. Dimeriza- 
tion is unconditional and increases Em iS) beyond E^ (0) 
as shown in Fig. [7] for a vertical excitation. The gs cusp 
appears again because the lowest triplet is not degener- 
ate and hence evolves as . The slope of (9£'m/9(5)o at 
the origin is NBiy^.N). It follows 



model is the generic case for e-ph coupling a — {dt/du)Q 
in a harmonic ID lattice with force constant k. Linear 
coupling is retained in generalizations [l^ that require 
numerical methods and include electron correlation, 
quantum fluctuations, nonadiabatic or 3D effects, spin 
chains and models with site energies. Quantum fluc- 
tuations are particularly important for small S that is 
easily reversed locally. Fluctuations reduce but do not 
wash out S ss 0.18 in the SSH model of polyacetylene [35| . 

Now Eq. [1] reads 

N 

H{d) = Hei+ ^ 5pt{a}p^^ap+i,„ + h.c) 

p— 1,(7 

+ (14) 
p 

with td = a^/fc and constant 5p = {—lyau/k in the gs. 
The distribution 5p is subject to the constraint 5p = Q 
for a fixed chain length. The gs energy per site is 

eriS) ^ eo{6) + 5y2ed. (15) 

A minimum at J ^ implies a dimerized gs. A global 
adiabatic approximation has 5{T) > for T < Tp, the 
Peierls temperature, while a local adiabatic approxima- 
tion leads to domains walls between regions of opposite 




"^•^0 0.02 0.04 0.06 0.08 O.I 

6 

FIG. 7: Ground state energy Eq in Eq. [15] and magnetic gap 
Em of an EHM with dimerization 5 in Eq. [2] and — 0.12. 
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that two domain walls that change the Sp pattern 
without changing l^ad to lower Ej^ , and reduced 

Em are readily found for iV = 12 or 16. Generalization 
of 2Ws{N) in Eq. [5] is more suitable for spin solitons 
between regions with opposite Seq. Such simulations are 
beyond the scope of the present study. 

An adiabatic approximation gives a dimerized gs with 
±Seq that has to be relaxed locally for solitons. X-ray 
detection requires that Seq not be too small compared to 
the zero point amplitude, (5^)^/^, of the Peierls phonon 
hiop in Eq. 1151 

(S^) = edhwp/2t (17) 

Typical values [s^ of fkjp lOOcm"^, ~ 0.30 and 
t « 1500cm~^ in organic stacks return ((5^)^/^ w 0.10 
at T = that increases with T. Finite £',„ and small 
Seq minimize contributions from thermal excitations, in 
contrast to SSH or correlated Peierls systems in which 
small 6eq necessarily implies small Em that vanishes at 
6 = 0. An arbitrarily small perturbation selects one of 
the degenerate gs of a BOW phase at K, but there is 
no long-range order in ID at finite T. Solitons lower the 
free energy and result in local adiabatic approximations. 
At low T, the mean separation i?(T) between solitons is 
large compared to their widths 2^. Each soliton can then 
be centered on any of R{T) sites. The soliton density in 
an extended system with R> 2^ and TV — >■ cx) is 

,(r).521.„,(_|) (18) 

where fc^ is the Boltzmann constant and W = Ws for 
spin solitons or Wc for charge solitons. The = 17 
results in Fig. 2] are not sufficient for estimating 2^, but 
indicate 2^ > 30 for spin solitons and 2^ « 30 for charge 
solitons, consistent with the SSH estimate [llj of 2^ « 15 
for larger 6 = 0.18. Since Ws < Wc, spin solitons are 
thermally accessible and the condition R{T) > 2^ holds 
up to p{T) w 1%. 

IV. SPIN SUSCEPTIBILITY 

The magnetic gap Em opens at the boundary Vg of 
the spin-fluid and BOW phases, and it does so very 
slowly at a Kosterlitz-Thouless transition [T]. DMRG 
with PBC provides an independent calculation [lOj 
showing that Em{V) opens at Vs = IMt for the EHM 
with U = At. The gap is only 0.023t a.t V = 2.0 before 
reaching 0.241t at V = 2.10 and 0.63t at V = 2.20, 
just in the CDW phase. We are interested in large 
Em near the boundary Vc of the BOW and CDW 
phases in FigU As noted above, Vc is a metallic 
point with bond orders p{V) in Eq. [5] close to 2/%, 
the band limit. The band limit returns Em = for 



N = 4n, when the degenerate orbitals at e = are half 
filled, and Em = Atsin{n/N) for N ^ An + 2. Open 
boundary conditions with tijy = have intermediate 
Em ■ The same pattern is seen in Fig. |S] for the EHM 
at y = 2.5, with smallest E„,{N) for N = 4n, tiN = 1 
and N = 4n-|-2, tu^ = —1. Quite unusually, Em 
increases with N at V = 2.5. DMRG calculations [l3| 
of Em with PBC have minimum Em{N) at Af « 30 
for V = 2.2. Exact Em{N) in Fig. [H] with tuq = 0, 
±1 and V = 2.2 decrease with N . Rapidly increasing 
Em at Vc is also seen [l3| for other potentials and has 
important implications for modeling the spin suscepti- 
bility. In particular, larger A^ is not automatically better. 

The molar magnetic susceptibility xm (T) allows quan- 
titative comparisons jl3i i25i . i37i] for organic ion-radical 
solids with small spin-orbit coupling and g-factors close 
to the free-electron value, g = 2.00236. We take Xm{T) 
to be the spin susceptibility after standard corrections for 
diamagnetism and impurities. The full spectrum of Eq. 
[T] is now required, and charge degrees of freedom vastly 
increase the number of states. We extend exact results 
for xm{T) up to A^ = 10. The partition function of the 
EHM in Eq. [T]with even A^ is 

N/2 

Qn{T) = 11(2^ + l)exp{-EsriN)/kBT) (19) 

S=0 r=l 

The singlet gs is the zero of energy, _Eoi(A^) = 0, and 
Esr{N) refers to the state r with spin S. The molar spin 
susceptibility is 

x{2S+l)exp{-Esr{N)/kBT) (20) 

where /.ts is the Bohr magneton and Na is Avogadro's 
number. Finite Em = En leads to xm(0) = 0, in 
contrast to finite XAf(O) for V < Vg in the spin-fluid 
phase. We set E^ = in the BOW phase and consider 
finite A" to be a coarse-grained approximation of an 
extended system with a dense spectrum for E > Em. 

Figured shows xm{T) up to ksT = t for A^ = 8 with 
PBC {tiM = 1) and N = 10 with tin = -l&tV = 2.10, 
2.15 and 2.20. Finite-size effects become negligible at 
high T, as has long been recognized in spin chains. In 
addition, the broad plateau around kpT t depends 
weakly on V , which simply refiects the narrowness of 
the BOW phase. By contrast, the T w behavior is 
governed by Em{V) and is very sensitive to V at constant 
U , t. Although t K, 10^ AT is a small electronic energy, 
experiment is limited to much lower T. The inset of 
Fig. |9] expands the relevant range. Finite Em suppresses 
Xm at low T after which xm is almost linear in T. The 
BOW phase of a frustrated spin chain "Ss"] without charge 
degrees of freedom also has an almost linear xm{T) and 
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FIG. 8: Magnetic gap Em of the A'^-site EHM, Eq. [U with FIG. 9: Temperature dependence of the molar spin suscepti- 

U = U and the indicated Vs. Boundary conditions with bility Xm{T) of iV-site EHM with [/ = 4t and V in Eq. [T]The 

t\N = and ±1 are discussed in the text. The DMRG energy curves in the inset are for the same U , V . 
gaps are from ref. 10], with minimum Em around A'^ ~ 30 
for V = 2.2t. 

structures [111 indicate regular TCNQ~ stacks. 

a broad maximum. 



A BOW phase with Em > has xm{T) w up to 
T PS Ti followed by a linear increase for T > Ti. For 
example, Rb-TCNQ(II) has flO^ Ti « 140 K that roughly 
fixes t{V) in the inset to Fig. [9] and completely specifies 
XmIT) of the EHM with U = 4i. Although xm{T) 
increases almost linearly to T = 300K, the magnitude 
at 300 K rules out the V > 2.10i curves in the inset 
while the V = 2.10 curve with /fcsT/t w 0.15 at 300 K 
fails for Ti. An EHM with J7 = 4t is not quantitative 
for Rb-TCNQ(H). Improved fits are possible for small 
6eq « 0.01 but such modeling also entails variation of 
Coulomb interactions. 

The EHM boundary Vs between the spin fluid and 
BOW phases has been difHcult to model and has been 
reported @ to be as high as Vs ~ 2.02t at U = it. 
Since XAf(O) is finite in the spin fluid phase and Em 
opens slowly for V > Vs, very different xm(T) are then 
calculated in the BOW phase [5|. The evolution of 
Em{V,N) in Fig. |H]with N and DMRG results support 
the original estimate |l| of Kj ~ 1.86< at U — At based on 
excited-state crossovers in Table ID Rapidly increasing 
Em{V) as V approaches the metallic point Vc is found 
in the BOW phase of related Hubbard models [lo| . 
Large Em is needed for the Rb-TCNQ(II) susceptibility 
as well as for K and Na-TCNQ at high T where X-ray 



V. DISCUSSION 

Rice jl^ recognized the possibility of measuring e-mv 
coupling constants gn from polarized IR spectra when 
Ci symmetry is broken on dimerization. Several groups 
[29! . Isoj extended the procedure to extracting transfer- 
able gn for selected 7r-donors and 7r-acceptors. K-TCNQ 
was a prime example j40l | of a crystal with dimerized 
TCNQ^ stacks at 300 K. Polarized mid-IR spectra 
show coupled ts modes that are shifted to the red from 
the corresponding Raman transitions [H, |4^. Powder 
Rb-TCNQ(II) has virtually identical IR transitions [41 1 
whose polarization along the stack has been confirmed in 
single crystals [1^, but strikingly different temperature 
dependence. Raman spectra and gn of TCNQ~ are 
expected to be the same, since solid-state perturbations 
are usually small [H, [s^] . 

The 100 and 295 K crystal structures of Rb-TCNQ(II) 
decisively indicate [T9I] a regular stack of TCNQ~ at in- 
version centers and interplanar separation R — 3.174A. 
The 295 K structure is in excellent agreement with 
previous data [43j . and the triclinic space group PJ is 
retained at lOOi^. Low R factors and examination of 
thermal ellipsoids at 100 K place a conservative limit on 
dimerization of R^ — i?_ < 0.05A [ll|. Yet negligibly 
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small Xm{T) below 140 K implies a large Em and IR 
data indicates broken electronic Ci symmetry. Broken 
Ci symmetry in a BOW phase accounts naturally for 
large £',„ near the CDW boundary and for IR intensity 
at K. We take up the different temperature dependence 
of the Rb salt. 

The IR intensity in Eq. [H] goes as {dP{V, A)/dA)l 
[2^[23| and increases with V in the BOW phase as shown 
in Fig. \5\ The intensities for a spin soliton in Fig. [5] are 
smaller and vanish at the center. Each soliton introduces 
a Ci center between regions with ±B(y) and reduces 
{dP{V,A)/dA)o over « 2^ sites. We approximate the 
temperature dependence as 



Iir{T) _ 1 



(21) 



Here ps{T) is the spin density given by Xm{T)/xc, 
where xc = NAQ^P^/'^kBT is the Curie susceptibility. 
The intensity is 50% lower at 2£^ps{T) — 1 where spin 
solitons overlap. Eq. [21] relates two measured quantities 
through 2^, as shown in Fig. [TU] for 2^ = 60, which 
is in the expected range. The Xm{T) data is for Rb- 
TCNQ(II) from ref. [ig and gives ps{T). The intensity 
ratio of the 722 cm~^ mode is from Fig. [5] of ref. [44| . 
with open and closed symbols on cooling and heating. 
Similar T dependence is seen for other mid-IR modes 
of crystals [42|. The fit supports a BOW phase inter- 
pretation and electronic solitons rather than a specific 
microscopic model or parameters. A microscopic model 
must account for xm{T) in addition to the intensity ratio. 

We have examined the BOW phase of the EHM at 
intermediate correlation [/ = 4i by direct solution of Eq. 
[T]for finite N. We used degenerate gs at y = Vi{N) in 
Table U to break inversion symmetry in finite systems 
and varied J2 in Eq. [5] to scan Vi{N) over the BOW 
phase. Exact degeneracy enforced by J2 turns out to be 
important for e-mv coupling to Holstein phonons but not 
for the order parameter B{V) = p+ — p- in Eq. [3l We 
compared BOW properties due to electronic correlation 
in a regular ID chain to the SSH model of a dimerized 
band. SSH results for topological solitons carry over for 
many aspects of spin and charge solitons in the BOW 
phase, albeit with different 6eff for different properties. 
We focused on the consequences of a degenerate gs 
and finite Em, especially large Em close to the metallic 
point Vc- Broadly similar results are expected in BOW 
phases of other Hubbard-type models with intermediate 
correlation. 

In the adiabatic approximation for the lattice, linear 
e-ph coupling generates a dimerized gs in both the 
BOW phase of the EHM and the SSH model, but they 
are different. The SSH model has a standard Peierls 
transition at Tp. Thermal population of excited states 
stabilizes the regular array for T > Tp, and low Tp 




100 150 200 250 300 



FIG. 10: Intensity ratio 1 1 r{T) / 1 1 r{0) of the 722 cm"^ 
mode of Rb-TCNQ(II) on cooling (open symbols) and heat- 
ing (closed symbols) from ref 44]. The solid line is Eg. 1211 
with 2^ — 60 and spin density ps{T) from Xm{T) in ref. [191 ] 



necessarily implies weak coupling or a stiff lattice. 
The Peierls transition of polyacetylene is far above its 
thermal stability, and we are not aware of evidence for 
5eq{T) variations up to ~ 400 K. Spin- Peierls systems 
illustrate decreasing 5eq{T) up to Tsp < 20 K that 
can be modeled in the adiabatic approximation {isj. 
The BOW phase of the EHM at intermediate U = At 
samples a different sector of parameter space, one in 
which substantial Em up to ~ 0.5t suppresses thermal 
excitations. Spin solitons in Eq. [T8| with Ws{S) give 
a small XAiiT) in this range, and Em = '^Ws remains 
finite at 5 — 0. The BOW phase has novel aspects that 
need further study. The principal theoretical issues are 
quantum fluctuations or nonadiabatic phonons that may 
suppress a sharp Peierls transition when e-ph coupling 
is weak. The experimental problem is to detect small 
dimerization against a background of zero-point motions. 

The present discussion is limited to the BOW phase 
of the EHM at U = At. We have developed the 
consequences of coupling to lattice phonons and to 
molecular vibrations in the adiabatic approximation. 
Similar results are expected [lo( for other quantum 
cell models with electron-hole symmetry and will be 
needed to model physical systems with BOW phases, 
starting with Rb-TCNQ(II). Since alkali-TCNQ sahs are 
semiconductors, they have Coulomb interactions rather 
than a Hubbard U and are stabilized close to the CDW 
boundary by the 3D electrostatic (Madelung) energy 



11 



[23 |. The Na and K-TCNQ salts have dimerization 
phase transitions [s^ with some 3D character since the 
cations also dimerize. The regular structure at high T 
has small xm{T) that increases with T in a manner that 
suggests a BOW phase. Both tt— radical organic stacks 
and conjugated polymers are quasi-lD systems whose 
initial modeling is without interchain interactions. 

In summary, we have characterized the BOW phase of 
the EHM with intermediate U = 4t by exact treatment 
of finite systems with degenerate gs at V = Vi{N). The 
elementary excitations are electronic solitons, both spin 
and charge, in a regular array. Solitons in the correlated 
BOW phase resemble the familiar solitons of the SSH 
model for e-ph coupling in a tight-binding band. Several 
measures indicate an "effective" dimerization 5 < 0.10 
that is considerably less than 5 = 0.18 for the SSH 
model of polyacetylene. Charge fluctuations are coupled 
to molecular (Holstein) phonons that become IR active 



in the BOW phase due to broken Ci symmetry. The T 
dependence of IR modes is consistent with spin solitons 
whose width is 2^ 60 lattice constants. The BOW 
phase is dimerized at K in the adiabatic approximation, 
but gs degeneracy and finite E„i lead to novel aspects 
for a possible Peierls transition. Previous discussions 
of SSH solitons greatly facilitate analysis of the BOW 
phase. So have previous treatments of e-ph and e-mv 
coupling in ID Hubbard models for conjugated polymers 
and organic ion-radical or charge-transfer crystals. The 
BOW phase of Hubbard-type models with intermediate 
correlation has some unique aspects that invite further 
study as well as features that are common to such models. 
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